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We propose the use of modulated spectra of astronomical sources due to gravita- 
tional lensing to probe Ellis wormholes. The modulation factor due to gravitational 
lensing by the Ellis wormhole is calculated. Within the geometrical optics approx- 
imation, the normal point mass lens and the Ellis wormhole are indistinguishable 
unless we know the source's unlensed luminosity. This degeneracy is resolved with 
the significant wave effect in the low frequency domain if we take the deviation from 
the geometrical optics into account. We can roughly estimate the upper bound for 
the number density of Ellis wormholes as n < 10 -9 AU -3 with throat radius a ~ 1cm 
from the existing femto-lensing analysis for compact objects. 
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I. INTRODUCTION 

In a variety of cosmological models based on fundamental theory, exotic astrophysical 
objects which have not been observed are often predicted. Conversely, an observational 
evidence for an exotic object would stimulate creative theoretical discussions. Probing these 
exotic objects and detecting them will give us significant progress of research in fundamental 
physics. Even if we cannot detect it, giving a constraint on the abundance of the exotic 
objects is one of powerful means of investigating the nature of our universe. 

Generally, the interaction between such unobserved exotic objects and well known matters 
is very weak or not well established. Thus, only the gravitational interaction would cause 
reliable observational phenomena. One of the most direct measurements of gravitational 
effects of an exotic object is gravitational lensing. For instance, massive compact halo 
objects are probed by using micro-lensing lJ-SL Cosmic strings are also targets for probing 
by using gravitational lensing phenomena 4Ml2|. In this paper, we propose a way to probe 
Ellis wormholes 13| by using lensed spectra of astronomical sources. 

The Ellis wormhole was first introduced by Ellis as a spherically symmetric solution of 
Einstein equations with a ghost massless scalar field. The dynamical stability of the Ellis 
wormhole is discussed in Ref. [l4j and the possible source to support the Ellis geometry 



was proposed in Ref. 15[. Gravitational lensing by the Ellis wormhole was studied in 



Refs. [16|, ll7[ and recently revisited by several authors [18|, ll9[. So far, it has been suggested 



that Ellis wormholes can be probed by using light curves of gamma-ray bursts 20|j, micro- 
lensing 2lM23| (see also Refs. 2J, |25[) and imaging observations 26|, |27j], while our proposal 
is the use of spectroscopic observations to probe Ellis wormholes. 

In order to fully investigate the lensed spectrum of a point source, the wave effect in 
gravitational lensing must be taken into account. The wave effect for the point mass lens 
is discussed in Refs. I28l . 1291 1. Wave effects in gravitational lensing by the rotating massive 

singular isothermal sphere 32J and the cosmic string 33|, [34 



object [3y], binary system 
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have been considered. In Sec. [TTl we calculate the amplification factor of gravitational lensing 
by the Ellis wormhole taking the wave effect into account. The geometrical optics limit is 
analytically presented in Sec. IIHI The difference in the amplification factor between the 
point mass lens and the Ellis wormhole lens is discussed in Sec. [IV] based on observables. 
In Sec. [Vj possible observations to probe Ellis wormholes are listed. Sec. [VI] is devoted to a 
summary. 

In this paper, we use the geometrized units in which the speed of light and Newton's 
gravitational constant are both unity. 



II. A DERIVATION OF THE LENSED WAVE FORM 



The line element in the Ellis wormhole spacetime can be written by the following isotropic 
form: 



(1) 
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where R = a corresponds to the throat and we simply call a the throat radius in this paper. 1 
Assuming the thin lens approximation is valid, we consider the wormhole lens system 
shown in Fig. Q3 We use the position vector X = (X, Y, Z) in the flat space. Then, the 

V s 




FIG. 1: Lens system with thin lens approximation. S, L, and O represent the source, lens, and 
observer positions, respectively. The path SAB is a ray trajectory which is specified with the vector 
£ on the lens plane £a- B' is the intersection of the line AO and the plane £3- £' is the position 
vector of the point B on the plane £b- 



coordinate R is given by R = \X — Xl|, where |Xl| is the lens position. We set Z-axis as 
the perpendicular direction to the lens plane and the source plane. -Ds, D^, and -Dls denote 
the distances from the observer plane to the source plane, from the observer plane to the 
lens plane, and from the lens plane to the source plane, respectively. 

In the geometrical optics limit, we consider light rays emanated from the source. The 
vector £ on the lens plane £a in Fig. [T] specifies the light ray which is deflected once at 
X = X L + £. Since £ := |£| can be regarded as the closest approach of the light ray, as is 
shown in Ref. 16(, the deflection angle a is given by 



■®"(i) ,+ °(i) 4 - 

As a result, the Einstein radius £o for the Ellis wormhole is given by 



(2) 



-7ra 2 D 



1/3 



(3) 



1 Since the throat surface area is given by 16ira 2 , our definition of the throat radius is half the areal radius 
of the throat. 
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where 



D 



(4) 



Since we are interested in the wave effect, that is, the deviation from the geometrical 
optics limit, we need to treat the wave equation rather than light rays. Neglecting the 
polarization effect, we consider the scalar wave equation with the frequency u. The wave 
equation for the monochromatic wave e lult <f)(X) is given by 



UJ 2 (j) ' 



2 \ -3 



,2 \ 



-AnA 5{X - X s ) 



(5) 



where X$ is the position vector of the point source and A in the source term is a constant 
which specifies the amplitude. Without the wormhole, we obtain the wave form 0o at the 
observer O as follows: 



.4 



cxp 



iuj\l -D| 



A 






iu)D s | 




(6) 



where t] = \rf\ and ff = (X s — X L , Y S — Y L , 0) and we consider the case 77 -C D s in this paper. 
Our assumptions to calculate the wave form at O are summarized as follows(see Ref. 29|): 



(a) The geometrical optics approximation is valid between the source plane and the plane 
E B in Fig. m 

(b) Thin lens approximation is valid and a ray from the source is deflected once on the 
lens plane Ea- 

(c) Assuming Ds ~ -Dl ~ -Dls ~ D, we use a non-dimensional parameter e defined by 
e := C,o/D, which gives the typical scale of the deflection angle. Then, we assume 
l/(uD) < e < 1 and rj/D = 0(e). 

(d) On the plane £3, the gravitational potential of the lens object is negligible and 
5D/D = 0(e), where SD is the distance between the planes Ea and Eb. 

On the assumptions made above, we calculate the wave form on the plane Eb up to the 
leading order for the amplitude and next leading terms for the phase part. Applying the 
Kirchhoff integral theorem between the plane Eb and the observer, we calculate the approx- 
imate wave form at O. 

The vector £ specifies the light ray which is deflected once on the plane Ea at X = X L + £ 
and reaches the plane Eb- The deflection angle is fixed by £ and the background geometry. 
The point A in Fig. [1] denotes the deflected point. We label the intersection of the deflected 
light ray and the plane Eb as B, while B' in the Fig. [1] denotes the intersection of the line AO 
and the plane Eb- As will be mentioned at the end of this section, the dominant contribution 
to the wave form at O comes from rays which satisfy £ ~ eD, where £ = |£|. Therefore we 
consider ^/D as 0{e) hereafter. 
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First, we consider the following ansatz for in the region between the source plane and 
the plane S B : 

<t> = f(X)e iS &. (7) 
On the plane E B , the amplitude f(X) is given by 

A 



[l + 0(e)). 



(8) 



In the geometrical optics approximation, the phase S(X) satisfies the eikonal equation given 
by 



FdiSdjS = uj z [1 + — 



,2 \ 2 



(9) 



At the point B, the phase based on the source position S is given by the following integral 



f dx i 
S| B = ^ -d.Sil, 

where we have introduced the optical path length / defined as 

dx l dxi 

Since, in the geometrical optics approximation, we find 



(10) 



djS = u[l 



R 2 



dx l 



(ii) 



(12) 



the integral ( fTUI) is given by 



S\ B = lu dl + cua 2 



R 2 



dl. 



(13) 



After the calculations explicitly shown in Appendix |AJ we finally obtain the following ex- 
pression: 



S\ B ^u 



D s 1 + 



2D 2 



+ 



2D LS \D L 



V_ 



ira 



(14) 



This expression for the phase and Eq. (jSJ) for the amplitude can be used for any value of £, 
that is, we have obtained an approximate wave form on the plane Sb- 

Applying the Kirchhoff integral theorem 35] and neglecting the contribution from the 
infinity, we express the wave form 0o at O by the following integral: 



1 

Air 



^B 



d /e 



'dZ V r 



e lujr d(f) B 
r dZ 



(15) 
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where 0b is the waveform at B. Since we are interested in only the leading order of the 
amplitude, we obtain 



iooA 



exp 



iuD* 1 



2DI 



d£ exp 



" w " 27rD L D LS 
where we have used the following approximations: 



ioj < -+ 



DlDs J_ V_ 
2D LS \D L D s 



ira 







dZ V r 



r 



IOJ 

A" 



8Z D LS 

Defining the amplification factor F by F := (fro/^o, we obtain 



(16) 

(17) 
(18) 



r ~ — / dx exp 

7ri 



where 



d 



x 



y 



e D S 



iud { (x — y) H — 



2^ 



2D LS D L 



(19) 

(20) 
(21) 
(22) 



c? gives the optical path difference between the lensed trajectory and the unlensed one in the 
geometrical optics limit for ff = 0. 

Introducing a polar coordinate, we rewrite this integral as 



„ ud 
F ~ — ( 

iri 



iujdy 2 



POO 


\- J 2 2\1 


/ dxx exp 


zojgs 1 i H — J 


Jo 





2- 



-2iude iu3dy2 



POO 




/ dxx exp 


jwd 1 i H — J 


/o 





d<p exp [— 2iuidxy cos </?] (23) 
J {2udxy). (24) 



The integrand is divergent at the infinity on the real axis. This is caused by our approx- 
imation associated with e, and not real. If we write down the integrand in a precise form 
without any approximation, we do not have any divergence. Actually, in the precise form, 
the contribution from the integral in the region x 3> 1 is negligible due to the cancellation of 
the quasi-periodic integration. For the same reason, the contribution from the integration in 
the region x <C 1 is also negligible. Since the approximate expression of Eq. (j24"j) is valid in 
the region i C 1/e, we can obtain a reliable result by neglecting the contribution from the 
region i > 1 in the integral f )24|) . Practically, to make this integral finite, it is convenient 
to consider the analytic continuation to the complex plane and take the path of the integral 
as shown in Fig. [2j Then, this integral can be numerically performed. 



7 




Rex 



FIG. 2: The path of the integral (|24p taken in the numerical integration. 



III. GEOMETRICAL OPTICS APPROXIMATION 



In this section, we derive an approximate form of the amplification factor F. In the 
expression ( l23l) . we apply the stationary phase approximation to the integral with respect 
to x. Then we obtain 



7T? 



2- 



f ~ \ l— rv 







v/r+2/ 



exp [zu;<i/i(so)] 



where 



/i(x) = x 2 — 2xy cos (p H — 

a; 



and s = ^o(^) > satisfies 



h'(x ) = <^ Xq — Xq?/ cos </? — 1 = 0. 



(25) 



(26) 



(27) 



Note that Eq. 0271) has only one positive root as a function of ip. 

In Eq. fl25|) . we again perform the stationary phase approximation in the integral with 
respect to <p, and we obtain 



F ~ F« 



-co 



V(4 + 2)(4-i) 



y^arl + 2)(1 -x 3 _) 



exp 



exp 



20; d 
iud 



-x\ + 4 

x + 

-xt +4 i7T 

si T 



where s± satisfies 



x| =F x|y -1 = 0. 



Note that 1 < x + and < s_ < 1. If we define /x± and #± as 

,,6 



(4 + 2)(4-i)' 



6L = ud 



—X± +4 7T _^ 7T 

si 4 4' 



(28) 
(29) 

(30) 
(31) 
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F geo can be expressed as 



goo 



(32) 



fi± is the magnification factor for each image in the geometrical optics approximation. 

As an observable, we focus on \F\ 2 in this paper. In the geometrical optics approximation, 
we obtain 

l-Fgeo| 2 = + + 2 v / | / u_ H u_| sm(2udr(y)) (33) 



with r(y) being the following: 



r(y) 



6_-9 + + vr/2 
2ujd ' 



(34) 



where note that the definition of r can be written in terms of x±, which is a function of y. 
\F\ 2 and l-Fgeol 2 are depicted as functions of u for each value of y in Fig. [3j 
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FIG. 3: \F\ 2 and |F gco | 2 for the wormhole lens, where T = dr. 



IV. COMPARISON WITH THE POINT MASS LENS BASED ON 

OBSERVABLES 

In this paper, we assume the following situations for the observation: 
• We can observe the spectrum of a source. 
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• The unlensed spectrum shape is well known. 

Note that, in our analysis, knowledge about the luminosity is not necessary. The amplifica- 
tion factor for the point mass lens is summarized in Appendix El For both point mass and 
wormhole cases, in the geometrical optics approximation, we obtain the form of Eq. ( |33l) or 
equivalently Eq. ( 1B3I) . 

In the frequency region where the geometrical optics approximation is valid, there are 
basically three observables which characterize the form of the amplification factor. The first 
is the frequency u, the second is the period of the oscillation of the spectrum as a function 
of oj, and the third is the ratio k between the amplitude of the oscillation and the mean 
value. The period of the oscillation of the spectrum makes T := r(y)d an observable, k is 
given by 



2y/\»+ 






+ 





(35) 



from Eq. f l33|) and plotted as a function of y as is shown in Fig. [H Since k is observable, y 



1 .u 



0.8 



0.6 




wormhole 
point mass 




0.0 



0.0 



0.5 



1.0 



1.5 

y 



2.0 



2.5 



3.0 



FIG. 4: k as a function of y for the Ellis wormhole case and the point mass lens case. 

and hence r(y) can be determined if we have enough accuracy of the observation. As shown 
in Fig. [5j r(y) is a monotonically increasing function of y and close to 2y in the region y < 1. 
Then, from another observable T = r(y)d we can obtain the value of d. 

The situation for the point mass lens case is the same as for the Ellis wormhole case. 
That is, the three observables u, T, and k can be regarded as gravitational lensing by a 
point mass as well as a wormhole. This fact indicates that we cannot distinguish which is the 
lens object only by using these three observables in the geometrical optics approximation. 
This degeneracy is resolved in the small frequency region in which the wave effect becomes 
significant as is explicitly shown in Fig. [6j 
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FIG. 5: r(y) for the Ellis wormhole case and the point mass lens case. The two cases are almost 
indistinguishable from each other in the region depicted here. 
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FIG. 6: Comparison between \F\ 2 and \F P \ 



To simply see this resolution of the degeneracy, we consider the small frequency limit, i.e., 
uj — > 0. In this limit, we have F — > 1 as explicitly shown in Figs. [3]and[8j If we can observe 
the spectrum in any frequency region of interest, the following quantity is an observable: 



lim IFI 2 



lim < \F\ > 

U)— >oo 



l//i tot := l/O+l + 



(36) 



where the bracket < > denotes the average through several periods. We depict 1/jitot as 
a function of y for both the wormhole case and the point mass lens case in Fig. [71 In both 
cases, 1/yUtot approaches to and 1 in the limits y — > and y — > oo, respectively. While 
V/^tot < 1 is satisfied in all domain of y for the point mass lens, l//xtot can exceed 1 for the 
wormhole case due to the demagnification effect originated from the negative mass density 
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FIG. 7: 1/fi 



tot- 



surrounding the Ellis wormhole. The behaviours of l//itot are totally different from each 
other. This fact shows that they are, in principle, distinguishable. 



V. OBSERVATIONAL CONSTRAINT 



One of the possible sources is gamma-ray bursts, which have been proposed to be used 
for probing small mass primordial black holes 36M39| and low tension cosmic strings 34 . 
Recently, the femto-lensing effects caused by compact objects were searched by using gamma- 
ray bursts with known redshifts detected by the Fermi Gamma-ray Burst Monitor 40J. From 
non-detection of the femto-lensing event, a constraint on the number density of compact 
objects have been obtained. For a fixed value of the mass M of each compact object, 
the constraint can be translated into a constraint on Qqo, where Qqo is the average energy 
density of compact objects in the unit of the critical density p cr . Then, the abundance of dark 
compact objects is constrained as Qqo < 0.15 at the 95% confidence level for M ~ 3 x 10 18 g. 

Since the number density n of the compact objects is given by 



n 



M 



2 x KT 9 AIT 



n 



CO 



M 



0.15 / \3x 10 18 g 



-l 



Pc 



10" 29 g/cm a 



(37) 



we obtain the constraint for the number dnsity of compact objects with M ~ 3 x 10 18 g as 
n<2x l(T 9 Ain 3 . Then, we can expect a similar constraint on the number density of Ellis 
wormholes with the throat radius which gives the same value of d as that for the compact 
object. The mass of 3 x 10 18 g gives d ~ 5 x 10~ 10 cm from Eq. (1B2I) and, from Eq. (|20|) . the 
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corresponding throat radius a is given by 



d 



) 



3/4 



( 



D 



) 



1/4 



a ~ 0.7cm 



5 x 10- 10 ' 



1Q28, 



(38) 



cm 



cm 



Therefore, the number density of Ellis wormholes with a ~ 1cm must satisfy n < 10~ 9 AU . 
Note that this constraint comes from the wave form in the geometrical optics approximation 
and hence we do not distinguish between the point mass lenses and Ellis wormholes. 

Another possible observation to probe Ellis wormholes is the observation of gravitational 
waves from compact object binaries. The unlensed wave form of the gravitational waves from 
a compact object binary is well known. From the chirp signal in the inspiral phase, we can 
obtain the spectrum of the gravitational waves. In order to distinguish the Ellis wormhole 
from a point mass lens, we need to observe not only the typical interference pattern but also 
the wave effect in the lensed spectrum. Hence, the spectrum in d ~ A := 2tt/uj is necessary 
to probe the Ellis wormhole. Assuming d ~ A, the typical throat radius of Ellis wormholes 
which can be probed by using gravitational waves is estimated as follows: 



The same estimate is applicable for galactic sources of electro-magnetic waves. We obtain 
a ~ 10 5 cm for galactic radio sources with A ~ 1cm, a ~ lm for galactic optical or infra-red 
sources and a ~ 1cm for galactic X-ray sources. The source must be compact enough to 
show the clear oscillation behaviour in the spectrum. This fact can be clearly understood 
by considering the y dependence of the phase in the amplification factor (13"3"|) . Since r(y) 
is roughly approximated by 2y, the period 5y for one cycle is given by 5y ~ n/(ud). The 
corresponding length scale 5r) on the source plane is given by 



where we have assumed d ~ A. If the source radius is larger than 5rj, the interference 
pattern will be smeared out. Observation of compact galactic sources such as pulsars and 
white dwarfs might be useful to probe not only dark compact objects but also exotic compact 
objects such as the Ellis wormhole. 



In this paper, we have proposed the probe of Ellis wormholes by using spectroscopic 
observations. We have assumed that the spectrum of the target source can be measured in 




(39) 




(40) 



VI. SUMMARY 
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enough accuracy and the spectrum shape is well known without lensing, but the luminosity 
is not necessarily observable. Then, we have discussed the distinguishability of the lensed 
spectrum from the case of the point mass lens. 

We have derived the wave form after the scattering by the Ellis wormhole including the 
wave effect in the low frequency domain. The geometrical optics limit of the wave form 
has been also analytically derived. Then, we have found that the Ellis wormhole cannot be 
distinguished from the point mass lens by using only the high frequency domain in which 
the geometrical optics approximation is valid. We have also found that this degeneracy is 
resolved in the low frequency domain in which the wave effect is significant. Possible obser- 
vational constraints are also discussed and we estimated the upper bound for the number 
density of Ellis wormholes as n < 1CT 9 AU -3 with throat radius a ~ 1cm from the existing 
femto-lensing analysis for compact objects. 

Finally, we note that our method to probe the Ellis wormhole is complementary to the 



other methods to probe the Ellis wormhole with micro-lensing 2lN23l] or the astrometric 



image centroid displacements [261] . These are not feasible for observations on cosmo logical 
scales because the time scale of the lens event is too long to detect modulation of the 
light curve or the displacements. In contrast, the slow relative motion is an advantage for 
spectroscopic observations. Therefore we may probe the Ellis wormhole on cosmological 
scales using our method. 
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Appendix A: Derivation of Eq. (JT 

The first term in Eq. (1131) can be evaluated as follows 

"B 



dl = ISAI + IABI = ISAI + I OA - O 



s 



|S^| + \/|0^| 2 + |OB>-2o1-OB\ (Al) 



OB 



and OB' can be written as 



oi 



, £ls + 6D f 
+ V+[£-V) „ ocSD- 



(l + 0(6 3 )) 



11 



0-f 








oi. 









4 



5D 



D 



LS 



;i + o(6 3 )), 



From these expressions, we can find 

|OB*| = |OB'| (1 + C(e 3 )) , 
OA* • Oil = OA ■ OB' (1 + C(e 3 )) 



Therefore we obtain 



Since we find 



|SA| = 

\A6\ = Dl (l 
we obtain the following expression: 



J dl = (\sA\ + \oA\ - |0F?|) (1 + C(e 3 )) 



(A2) 
(A3) 

(A4) 
(A5) 

(A6) 



2^8 



2^ 



(A7) 
(A8) 



£>s 1 + 



2D£J 2D LS \D L D s 



(1 + (9(e 3 )) , (A9) 



where r = |0B'|. 

In order to evaluate the second term in Eq. (I13j) . we first consider the integral between S 
and A. Letting P be a point on the segment SA, we obtain 



|LPT = || + AT^I 2 = e + ( 1 - -M Al 



e + (Isa -l) 2 + 2(l--^-) £ • AS, (A10) 



where Zsa = |SA| and / = |SP|. Since £-AS = —£•(£ — 77), we obtain the following expression: 

0\ 2 = e + a SA - n 2 - 2 (1 - J-] £■ d- ft). (ah) 



e 2 + (/ S A-o -2(^i-^je-(e-^). 

Substituting the above expression of \l3\ 2 into R 2 of the second integral in Eq. (fi~3l) with 
the integral region being from S to A, we obtain 



SA 



dl 



e + (isA-iy-2(i-^-(t-ff) 



(A12) 
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Then the integral flA12l) can be performed and evaluated as 



is A 



dl 



(lsA-l) 2 + e-2(l-^-(^-v) 6Dls 



arctan 



7T £-f) 



2£ ?D 



LS 



-*SA + Z-(Z-fl)+ USA 

(l + 0{e 2 )). (A13) 



The contribution from the integral between A and B can be also evaluated by the similar 
integral. Finally, we obtain the expression ( |T4|) . 



Appendix B: Point Mass Lens 
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For the point mass case, we have the following expression for the amplification factor 28 



| F po| = | e ™d/2 r ^ _ iu(£) ^ ( iud ^ j. iudy j\ 



(Bl) 



where T and \F\ are the gamma function and the confluent hyper-geometric function, re- 
spectively, and 



d = 2M, y = r] 



AMD LS D S 

In the geometrical optics approximation (u — > oo), we obtain 



where 



I pp°| 2 _s. I pp° I 2 

I gco 



po 



1/4° I + Im— °l + 2 V l/4V-°l sin(2o;drp (j/)), 



1 

±- 

4 



± 2 



T po(j/) 



vV + 4 + y 



2 ^2 + 4 _ y 



(B2) 
(B3) 

(B4) 
(B5) 



F po | 2 and |i^ c °| 2 are depicted as functions of u for each value of y in Fig. |HJ 
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